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Abstract 

Modeling of headway /spacing between two consecutive vehicles has many applications in traffic 
flow theory and transport practice. Most known approaches only study the vehicles running on 
freeways. In this paper, we propose a model to explain the spacing distribution of queuing vehicles 
in front of a signalized junction based on random-matrix theory. We show that the recently mea- 
sured spacing distribution data well fit the spacing distribution of a Gaussian symplectic ensemble 
(GSE). These results are also compared with the spacing distribution observed for car parking 
problem. Why vehicle-stationary-queuing and vehicle-parking have different spacing distributions 
(GSE vs GUE) seems to lie in the difference of driving patterns. 
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I. INTRODUCTION 



The spacing is usually denned as the distance between two successive vehicles measured 
from the same common feature of the vehicles (e.g. rear axle, front bumper). Because the 
distribution of spacing reflects the unmeasurable interaction forces or potentials between 
vehicles that governs their motions, increasing investigation are put into this field to reveal 
the complex dynamics of vehicle traffic flow and explain some important phenomena, i.e. 
phrase transition {2}, Q], 0]. 

One interesting topic is to discuss spacing distributions observed during the formation of 
and transitions between different vehicle queues: static queues (vehicles parked on a line in 
the parking lot, or vehicle queues fully-stopped in front of signalized intersections), moving 
queues which may contain diversified inter-arrival and inter-departure queuing interactions 



0, Ei. 



In this short paper, we will focus on the not so popu 



ar static queues. 
7] as: how many randomly 



The vehicle parking problem was first introduced by Renyi in 
parking motorists can be accommodated on a line street of a given length on average. The 
most famous solution to this question is based on Random Sequential Adsorption (RSA), 
which looks it as an irreversible process in which particles are adsorbed sequentially and 
without overlap onto randomly chosen positions on a surface. This ID RSA problem can 
be solved analytically when all the vehicles with the same equal length, see 
However, it is hard to directly apply this method to other vehicle queues by definition. 

Differently in H[, the random-matrix theory is used to study the car-parking problem, 
where the nature of interaction between the particles in a Dyson's Coulomb gas model is 
assumed to be consistent with the tendency of the drivers to park their cars near to each 
other and in the same time keep a distance sufficient for departure. It was shown that 
the recently measured gap-size distribution of parked cars in a number of roads in central 
London can be well represented by the spacing distribution of a Gaussian unitary ensemble. 

For a similar purpose, in this paper, the Dyson's Coulomb gas model is adopted to explain 
the formation of vehicle queues fully-stopped in front of signalized intersections. We found 
a good agreement between the empirical data and the spacing distribution for Gaussian 
symplectic ensemble of random matrices. 
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II. COULOMB GAS MODEL AND QUEUING VEHICLES 



Considers a gas of iV charges whose positions are denoted by x±, £2, %n and these 
charges are free to move on the line < x < +00, see Fig. [TJ Suppose the potential energy 
of this Coulomb gas is given as [12] 



V 



(1) 



i<j 



where the first term in (TjQ) represents a harmonic potential that attracts each charge 
independently towards the coordinate origin; and the second term represents an electrostatic 
repulsion between each pair of these charges. 
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FIG. 1: Queuing vehicles in front of a signalized junction analogue to the Dyson Gas. 

Let P(s) denote the nearest-neighbor spacing distribution of these charges. The accurate 
solution of P(s) is not easy to find. However, for such systems, the probability density 
function for the position of the charges can be approximately calculated by using the so 
called Wigner surmise [3], [l4], 15 1. 



Suppose the gas is in thermodynamical equilibrium at temperature T = 1/(K(3), where 
K is the Boltzmann constant. The probability density function can be given as below by the 
Boltzmann factor obtained from the Gibbs-Boltzmann canonical distribution by integration 
over the momenta of the particles. 



P(xi,x 2 , -:,x N ) = Ce 



-pv 



(2) 
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Combining ([I]) and (J2J), the Wigner surmise solutions for P(s) can be gotten by taking 
the additional assumption of (3. The role of this inverse temperature f3 denotes the level- 
repulsion power of the matrices eigenvalues. Particularly, they are suggested by Wigner as 
below respectively. For (3 = 0, we get the well-known Poisson Ensembles (PE) 

Ppe(s) = e~ s (3) 
for (3 — 1, we get Gaussian Orthogonal Ensembles (GOE) 

Pgoe{s) = ivrse-^ 2 (4) 
for (3 = 2, we get Gaussian Unitary Ensemble (GUE) 
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Pgue(s) = -^e~^ (5) 
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and for f3 = 4, we get Gaussian Orthogonal Ensembles (GSE) 
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P ose(s) = ^e'^ (6) 

Similar as discussed in [111 ], we can abstract the movements of vehicles of different size 
into point particles, because we are only interested in the spacing distribution here. A 
natural guess for the queuing dynamics of vehicles in front of a signalized junction is that 
the system can also be approximately formulated into this Coulomb model. 

The single-particle term in ([1]) can be viewed to reflect the tendency of driving closer and 
the repulsive two-body term in (CQ) indicate the tendency of maintaining the safe distance. 
Analogously, the basic instinct of a driver is to maintain a small and safe gap between 
him/her and his/her leading vehicle, especially when he/she is queuing. The superposition 
of these two potentials, which appears an overall repulsion for small spacings and attraction 
for the large ones, expresses the fact that it is unlikely to see too small or too large spacings 
between queuing vehicles. However, no one can always keep an ideal headway due to distur- 
bances (unexpected acceleration/deceleration of the leading vehicle, occasional absence of 
mind, etc.). Thus, these vehicles (particles) are perturbed by environment simultaneously. 

Noticing the above analogue and inspiriting by the report that the spacing distribution of 
vehicle parking is in agreement with GUE type distribution, we conjecture that the empirical 
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spacing distribution of a queuing vehicle systems might fit one kind of Wigner surmises (J3J)- 
©, too. 

III. COMPARISON WITH EMPIRICAL RESULTS 

To test this conjecture, we collected 700 sample spacings of queuing vehicles in front of 
several different signalized junctions in Beijing, China. Some details about data collection 
can be found in [161 ]. 

The average spacing size observed is 1.43m here. Fig. [2] shows the probability distribution 
function P(s) in a form of normalized histogram, where the values of the x-axis is defined 
as the ratio of the spacing to the mean value. This modification is introduced to make the 
data comparable with the nearest-neighbor spacing distribution for system (CQ). 




Data 

FIG. 2: The normalized spacing distribution of queuing vehicles in front of several signalized 
junctions in Beijing compared with the theoretical spacing distribution for Poisson, GOE, GUE, 
GSE. 

In the Coulomb gas model, the inverse temperature f3 of the gas characterizes the degree 
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of repulsion. Fig. [2] shows the theoretical spacing distribution curves for Ppe(s), Pgoe(s), 
Pgue(s) and Pgse(s) as well as the empirical spacing distribution histogram. We can see 
that different from the vehicle-parking scenarios, the vehicle- queuing scenarios well fit the 
GSE type model instead of GUE type model, although in these two scenarios, drivers all aim 
to driving close enough but not too close. This suggests that the vehicles queuing process at 
a signalized junction can be added to the long list of the systems with RMT-like fluctuations. 

An interesting question is why the spacing distributions of vehicle-parking and vehicle- 
queuing are different. A reasonable explanation is this difference comes from the dissimilar 
driving patterns. In vehicle-parking scenarios, drivers would like to try several times and 
move back-and-forth to adjust the gaps so as to park to an "ideal" position; while in vehicle- 
queuing scenarios, driver won't be able to drive back. This phenomenon can be interpreted 
as: in the vehicle-parking scenarios, the repelling force from the neighboring vehicles is 
relatively "loose" ; however, in the vehicle-queuing scenarios, the repel force from the neigh- 
boring vehicles is a kind of "rigid" . As pointed out in [7] , such a difference will result in the 
different values of inverse temperature (3. At low temperatures ( (3 is bigger), the charges 
tend to be regularly spaced in a crystalline lattice arrangement and the randomness of the 
positions of the charges are small. At higher temperatures ( (3 is smaller), the fluctuations 
of the charges become tenser. Thus, for vehicle-parking scenarios, we get (3 = 2 and for 
vehicle-queuing, we get (3 = 4. 

We also guess that in different cities, the spacing of vehicle queues may still hold GSE 
type distribution but may have different mean value. Further experiments will be carried 
out to test this guess soon. Any vehicles queuing data collected in cities other than Beijing 
are welcome and highly appreciated. 

Acknowledgments 

This work was supported in part by National Basic Research Program of China (973 
Project) 2006CB705506, Hi- Tech Research and Development Program of China (863 
Project) 2006AA11Z229, and National Natural Science Foundation of China 60374059, 



6 



50708055. 



[1] D. Chowdhury, L. Santen, A. Schadschneider, Phys. Rep. 329, 199 (2000). 

[2] D. Helbing, Rev. Mod. Phys. 73, 1067 (2001). 

[3] R. Mahnke, J. Kaupuzs, I. Lubashevsky, Phys. Rep. 408, 1 (2005). 

[4] M. Schonhof, D. Helbing, Transp. Sci. 41, 135 (2007). 

[5] G. Zhang, Y. Wang, H. Wei, Y. Chen, in Transportation Research Board Annual Meeting CD, 
2007. 

[6] D. Helbing, M. Treiber, A. Resting, Physica A 36, 62 (2006). 

[7] A. Renyi, Sele. Trans. Math. Stat. Probab. 4, 203 (1963). 

[8] J. Talbot, G. Tarjus, P. R. Van Tassel, P. Viot, Colloids and Surf. A 165, 287 (2000). 

[9] J. W. Lee, Physica A 331, 531 (2004). 

[10] S. Rawal, G. J. Rodgers, Physica A 346, 621 (2005). 

[11] AY. Abul-Magd, Physica A 368, 536 (2006). 

[12] F. J. Dyson, J. Math. Phys. 3, 140 (1962a). 

[13] T. Guhr, A. Muller-Groeling, H. A. Weidenmuller, Phys. Rep. 299, 189 (1998). 

[14] M. L. Mehta, Random Matrices, 3rd edition, (Academic Press, Boston, 2004). 



[15] M. Dieng, C. A. Tracy [arXiv:rnath/0 603543yl 



[16] Y. Su, Z. Wei, S. Cheng, D. Yao, Y. Zhang, L. Li, Z. Zhang, Z. Li, in Transportation Research 
Board Annual Meeting CD, (2008). 



7 



